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Abstract 



Standard derivations of "time-independent perturbation theory" of quan- 
CN ■ turn mechanics cannot be applied to the general case where potentials are 

', energy dependent or where the inverse free Green function is a non-linear 

| function of energy. Such derivations cannot be used, for example, in the con- 

text of relativistic quantum field theory. Here we solve this problem by pro- 
■ viding a new, general formulation of perturbation theory for calculating the 

changes in the energy spectrum and wave function of bound states and reso- 
nances induced by perturbations to the Hamiltonian. Although our derivation 
is valid for energy-dependent potentials and is not restricted to inverse free 
O j' Green functions that are linear in the energy, the expressions obtained for the 

energy and wave function corrections are compact, practical, and maximally 
similar to the ones of quantum mechanics. For the case of relativistic quan- 
tum field theory, our approach provides a direct covariant way of obtaining 
corrections to bound and resonance state masses, as well as to wave functions 
that are not in the centre of mass frame. 

ll.10.St, 11.80.Fv, 13.40.Ks, 31.15.Md 



Typeset using REVTgX 



*On leave from Mathematical Institute of Georgian Academy of Sciences, Tbilisi, Georgia. 

1 



I. INTRODUCTION 



There is a growing interest in calculations, within a covariant quantum field theory 
framework, of changes in the properties of bound states and resonances induced by small 
perturbations in the interaction Hamiltonian. The four-dimensional Bethe-Salpeter equation 
and its various three-dimensional reductions (so-called quasi-potential equations) are the 
most popular tools in this respect. A current example is the Nambu Jona-Lasinio (NJL) 
model where the nucleon is described in terms of three relativistic quarks interacting via 
contact potentials, and where meson exchange provides an important perturbative correction 
[1]. Another example is provided by relativistic calculations of hadronic atoms where the 
strong interaction perturbs the Coulomb bound state [2,3], and yet another by various other 
corrections to relativistic calculations of electromagnetic bound states [4]. 

The perturbation problem involved in such covariant calculations can be formulated as 
follows. Denoting the total four-momentum of the system by P, one would like to determine 
the bound state solution of the equation 

^ = (1) 

where K\(P) is a perturbation to the unperturbed kernel Kq(P), and where it is assumed 
that the unperturbed Green function G U (P), defined as the solution to the equation 

G U (P) = G (P) + G (P)K (P)G U (P), (2) 

is known completely. 1 Thus we seek the mass M and wave function \1/ such that Eq. (1) with 
P 2 = M 2 is satisfied. A consequence of the complete knowledge of G U (P) is that the mass 
spectrum M™ {n = 1, 2, 3, . . .) and corresponding wave functions of the unperturbed 
equation 

G \P)-K (P) 

where P 2 = (M™) 2 , are known. 

The task of solving Eq. (1) by expressing the mass M and wave function \l> as a perturba- 
tion series with respect to K\ is a problem whose solution is well-known in the corresponding 
context of non-relativistic quantum mechanics (given by so-called time-independent pertur- 
bation theory). Unfortunately the (textbook) derivation used to obtain the quantum me- 
chanical result is restricted to the case where the inverse free Green function Gq 1 (P) is lin- 
early dependent on energy Pq and where the unperturbed kernel Kq is an energy-independent 
Hermitian operator. Although these restrictions lead to the closure and orthonormality con- 
ditions 

$n$m = S nm , ®n$n = 1, (4) 

n 



For simplicity of presentation we generally do not show spin or relative momentum variables; 
similarly, identical particle factors and all sums and integrals over intermediate state variables are 
suppressed. 
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which are crucial for the derivation of time-independent perturbation theory, they are not 
valid in the Bethe-Salpeter case. Indeed none of these restrictions are required in the context 
of a covariant field theoretic approach. In this paper we therefore present a new solution 
to the perturbation problem which is valid for any form of Gq 1 (P) and Kq(P); in particu- 
lar, our solution is valid for the case of covariant field theoretic approaches where Gq 1 (P) 
depends nonlinearly on Pq and where Kq(P) can be energy (Pq) dependent. Our solution, 
given in Eq. (24) and Eq. (25) for the nondegenerate case, and in Eq. (43) and Eq. (46) 
for the degenerate case, expresses the mass M of the bound state or resonance and the 
corresponding wave function \1/ in terms of compact expressions that take into account the 
perturbation term K 1 to any order. At the same time, our formulation allows us to write 
the perturbation series for both M and up to any order, in a straightforward way which 
is maximally close to the analogous quantum mechanical formulation. A further important 
aspect of our approach is that it is manifestly covariant. This feature enables the direct use 
of the perturbation series for \1> also in cases where the bound state or resonance is not at 
rest. In this way the more involved approach of Lorentz boosting wave functions calculated 
perturbatively in the rest frame, can be avoided. As such, our approach to the perturbation 
problem where no restriction is put on the energy dependence of kernels and inverse free 
Green functions, may provide some important advantages over previous formulations [5,6,2]. 



II. PERTURBATION THEORY 
A. Basic equations 

In this paper we use the framework of relativistic quantum field theory to illustrate our 
approach to perturbation theory. Although this is done partly for presentational purposes - 
it is a particular case where the kernel is energy dependent and where the inverse Green func- 
tion is non-linearly dependent on energy, it is also a particularly topical discussed 
in the Introduction. On the other hand, we emphasize that our approach to perturbation 
theory does not depend on the particular theoretical framework in which the bound state 
problem is set - it can be that of non-relativistic quantum mechanics, relativistic quantum 
field theory, three-dimensional relativistic quasi-potential equations, etc. Similarly, our ap- 
proach does not depend on the functional form taken by the energy dependence of either 
the kernel or the inverse free Green function. All we need to assume is the usual overall 
structure of the dynamical equations involved, as exemplified by Eq. (1) and Eq. (2). 

We thus consider the Green function 

G(P) = G (P) + G (P)K(P)G(P), (5) 

where P is the total four- momentum, Go is the fully disconnected part of G, and where the 
kernel K consists of a part Kq for which the corresponding Green function is known, and a 
small part K\ which can be treated as a perturbation. Thus 

K(P) = K (P) + K 1 (P), (6) 

and it is assumed that the unperturbed Green function G U (P) has been previously deter- 
mined by solving Eq. (2). We are interested in the case where G U (P) has a pole corresponding 
to a bound or resonance state. Thus we can write 
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where the wave functions $(P) and $(P), the unperturbed bound state mass M u , and the 
background term G b u (P) are all assumed to be known. 2 In this respect it is worth noting 
that the pole term of Eq. (7) is separable with respect to initial and final state variables, 
thus for a two-body system $(P) = 3>(P,p) is a function of the initial relative momentum 
p while $(P) = <&(P,p') is a function of the final relative momentum p'. Note also, that 
as P — > P u , where P u is any four- vector such that P 2 = M 2 , the wave functions $(P) and 
$(P) must reduce to the respective solutions of the bound state equations 

$(P„) = G (P U )K (P U )<S> (P u ) ; $(P U ) = $(P U )K (P U )G (P U ). (8) 

Although $(P U ) and $(P U ) are therefore specified as the solutions of the above bound state 
equations, for momenta P not on the mass shell, P 2 ^ Ml, there is no unique way to 
define $(P) [and therefore $(P)] since any definition can be adopted in Eq. (7) with an 
appropriate redefinition of the background term G b u (P). Here we shall choose $(P) to be a 
Lorentz covariant function of the total momentum P, the relative momenta, and the spinor 
indices of the constituents (i.e. $(P) is covariant under the simultaneous transformation of 
all these variables). The way to construct such a $(P) will be discussed below. Since the 
full unperturbed Green function G U (P) is a Lorentz covariant function of its variables from 
the outset, the Lorentz covariance of the background term G b u (P) is therefore assured. 

Once the perturbation K\ is included, the mass M u will shift to the physical value M 
and $(P) will modify to the wave function \&(P) where 

G{P) jmm + &{P) . (9) 

The wave functions ^(P) and ^(P) are likewise assumed to be covariant functions which 
reduce in the limit P — > P, where P 2 = M 2 , to the respective solutions of the bound state 
equations 

tf(p) = G (P)K(P)V(P), and *(P) = V(P)K(P)G (P). (10) 

To write a perturbation series for C, we express G in terms of the known unperturbed 
Green function G u through the equation 

G(P) = G U (P) + G U (P)K 1 (P)G(P), (11) 

which follows from the fact that G^ 1 = Gq 1 — K and G^ 1 = Gq 1 — K . By iterating Eq. (11) 
we obtain a perturbation series for G(P) with respect to the perturbation K\{P). What 



2 Here, for simplicity, we assume that the bound state is nondegenerate - the degenerate case will 
be considered in detail in the next subsection. Also, here and elsewhere, all references to a "bound 
state" should be understood to include the case of a "resonance state" . 
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appears more difficult is to find a corresponding perturbation series for the mass M and 
wave function Yet if one closely examines the structure of the above equations, it can 
be discovered that a mathematically similar problem was solved long ago by Feshbach [7] 
albeit in the rather different context of nuclear reaction theory. Indeed there are a number 
of other contexts where analogous problems have been solved, the case of mass and vertex 
renormalization in pion-nucleon scattering being particularly noteworthy [8]. In the next 
section we shall therefore use the method of Feshbach to derive the solution of our covariant 
perturbation theory problem. 

B. Solution 

In this subsection we derive expressions for the bound state wave functions and the 

bound state mass M corresponding to the full kernel K of Eq. (6). Although our goal is to 
formulate the covariant perturbation theory for this problem, we in fact derive expressions 
for ty, and M, that are exact with all orders of K\ being taken into account. Starting 
from these exact expressions it is then trivial to generate all terms of the perturbation 
series. To present our solution it will be convenient to discuss the cases of nondegenerate 
and degenerate states, separately. 

1. Nondegenerate case 

In the nondegenerate case, to each unperturbed bound state mass M u there corresponds 
a unique bound state wave function $. The unperturbed Green function G U (P) then has the 
"pole plus background" structure, as given in Eq. (7). Having in mind that the full Green 
function G(P) has a similar structure as given in Eq. (9), and that our goal is to relate the 
quantities in these two expression, we begin by introducing a "background" Green function 
G b (P) defined as the solution of the equation 

G\P) = G b u {P) + Gt(P)^(P)G 6 (P). (12) 

Note that G\P) ^ G b (P) where G b {P) was defined in Eq. (9). From Eq. (12) it follows 
that 

(l + G b K 1 )~ 1 G b = G b u , (13) 

where we have dropped the momentum arguments for convenience. Similarly Eq. (11) implies 

G{\ + K 1 G)~ 1 = G n . (14) 
Subtracting the last two equations, we obtain 

G(l + KiGY 1 - (1 + G b K 1 )~ 1 G b = (15) 

and therefore 
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1 + G b K 1 )G - G\l + K ± G) = (1 + G b K x ] 



p2 _ M 2 



(l + K.G). 



Thus 



G = G + 



[1 + + KiG) 

P 2 - Ml 



which can be solved for $(1 + K\G) by writing 
$(l + ^iG) = $(l + ^iG 6 ) + 

and then 



,,, t <5>K 1 (l + G b K 1 )i<5><5>(l + K r G) 
P 2 - Ml ' 



§{\ + K x G) 



1 - 



P 2 -M 2 

Using this result in Eq. (17) we obtain the result we are seeking: 

i^(P)i>(P) 



G(P) 



P 2 - Ml - i$(P) [K^P) + K l (P)G b (P)K 1 (P)} $(P) 
where the functions ip(P) and ip{P) are defined by 

V>(P) = [l + G b (P)K 1 (P)] $(P) 

and 

^(P) = <5(P) \l + K 1 (P)G b (P) 



+ G b (P), 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



A comparison of Eq. (20) with Eq. (9) shows that ^(P) = y/Zi/>(P), and V(P) = ^/Zip(P), 
where 



Z = 



1 



(23) 



p2=p2 =M 2 



1 - i {<5(P) [K^P) + K 1 (P)G 6 (P)^ 1 (P)] $(P)}' 

with the prime indicating a derivative with respect to P 2 , and 

M 2 = Ml + il>(P) [Xi(P) + Ki(P)G 6 (P)ir 1 (P)] $(P). (24) 

In this respect it is worth noting that because all our wave functions and Green functions 
are Lorentz covariant, the quantity in the curly brackets of Eq. (23) [which also appears in 
Eq. (24)], is a Lorentz scalar depending only on P 2 . 

Thus, in the nondegenerate case, the properly normalized wave functions for the full 
perturbation theory are 

*(P) = {l - i {HP) [MP) + K 1 (P)G b (P)K 1 (P)] $(P)}'} V2 [1 + G 6 (P)K 1 (P)]$(P), (25) 



q?(p) = $(P)[1 + K 1 (P)G b (P)} |l - ^ {$(P) [k^P) + K l {P)G b {P)K l {P)] $(P)} 

We note that these wave functions satisfy the normalization condition 

dG~\P) 



-1/2 



(26) 



itf(P)- 



<9P 2 



-¥(P) 



= 1. 



(27) 



p=p 
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2. Reference frame dependence of the wave functions 

As far as we know, all previous attempts at developing perturbation theory for relativis- 
tic systems have considered bound states only at rest (see e.g. [5]). On the other hand, 
for observables involving scattering off the bound state (e.g. electromagnetic form-factors) 
taking into account the total momentum dependence of the bound state wave function is 
important. In the relativistic case there are some subtleties in the determination of this de- 
pendence perturbatively and at the same time in a manifestly covariant way. One possible 
way to do this is to derive the wave function to the needed order in the rest reference frame, 
and then to boost it in order to give it the desired momentum. There are two disadvantages 
to this approach: one is that it involves two separate steps - the perturbation expansion 
and the boosting. The second disadvantage is that the unit vector n = P/M which deter- 
mines the boost [14], itself may need to be calculated perturbatively. To illustrate this, we 
consider the determination of a scalar bound state wave function ^f(P) to first order in the 
perturbation. Showing explicitly one relative momentum p in addition to the total on-shell 
momentum P, we first write the perturbed wave function as a boosted wave function at rest: 

*(P,p) = S Ln V(L n P,L nP ) = S Ln V (L n p) (28) 

where L n is the boost Lorentz transformation, L n P = (M, 0), Si n is the associated trans- 
formation matrix acting on the spin indices of the constituents, and ^o(q) is the bound 
state wave function at rest. Next step is to calculate ^o(q) to first order in the perturba- 
tion: ^fo(q) = (1 +77i)$ (g), where the first-order correction factor rji is given explicitly in 
Eq. (57). Thus 

^(P,p) = S Ln {l + ri 1 )ML n p). (29) 

As L n is a function of the unit vector n = P/M = (VP 2 + M 2 , P)/M, and therefore of M, 
and because we need ty(P,p) up to first order, the mass M, should be approximated up to 
first order in the perturbation. Denoting the first order perturbation correction to M 2 by 
Si [given explicitly in Eq. (53)], the approximation n(M) xs n(M u + 8\/2M u ) should thus be 
used in Eq. (29) with a subsequent expansion of the resulting ty(P,p) up to first order in 5±. 
If admixtures of higher-order corrections were acceptable, then this last expansion could be 
neglected. 

In what follows we show a more straightforward way to obtain the perturbed wave 
function ^(P) when P^O. For this purpose we shall require G b (P), which determines the 
wave function via Eq. (25), to be Lorentz covariant; that is, we would like it to transform 
kinematically under any Lorentz transformation L of the momenta involved, as G b (P; p', p) = 
SLG b (LP; Lp', Lp)S\ j where p and p' are the initial and final relative momenta. In order for 
this to be satisfied, G b u {P\p' ,p) should also be Lorentz covariant in view of Eq. (12). Using 
the definition (7) for G b u (P;p' ,p) one can see that the unperturbed wave function $(P, p) 
should also be a Lorentz covariant function under any Lorentz transformation L of P and 
p. 

Thus the essential problem is a practical one: how to construct a wave function $(P,p) 
that is Lorentz covariant, and which satisfies the bound state equation [first of Eqs. (8)] for 
any P such that P 2 = M 2 . For this purpose it is useful to have a separate notation for the 
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bound state wave functions, so to this end we denote by $(P u ,p) all the solutions of the 
bound state equation [first of Eqs. (8)] for which the total momentum P u has the property 
Pj = M 2 . We then note that one cannot simply define <3>(P, p) = §(P u ,p) where P = (Po, P) 
and P u = (\JP 2 + M£, P), so that $(P,p) does not depend on Po - such a $(P,p) cannot be 
Lorentz covariant since a Lorentz transformation will change this function to Sl&{LP, Lp) 
which will necessarily depend on the (arbitrary) value of Po (the three-vector part of LP 
depends on P ). 

To make progress, we note that the bound state wave function $(P u ,p) is covariant under 
the transformation P u — > LP U , p — > Lp: 

$(P u ,p) = S L $(LP u ,Lp). (30) 

As this is true for any four- vector P u satisfying P 2 = M„, it will certainly be true for the 
four- vector M u P/\fP* where P is arbitrary. Thus, if we define wave function $(P,p) as 

*(P,p) = $^,p), (31) 

it immediately follows that 

^(P,p) = S L ^(LP,Lp), (32) 

which is the statement that wave function $(P, p) is Lorentz covariant in the way we need. In 
this way we have constructed a wave function $(P, p) that satisfies the sought-after Lorentz 
covariance, while at the same time reducing to the bound state wave function $(P u ,p) as 
P — > P u (in fact $(P, p), as defined by Eq. (31), is the bound state wave function with total 
momentum M U P/ y/P^). By choosing the form of $(P) given in Eq. (31), we guarantee that 
Eq. (7) is expressed in a manifestly covariant way. The immediate consequence of this is 
that the exact wave function ^(P) is given, up to a scalar normalization, in a manifestly 
covariant way by Eq. (21), and so is each term in Eq. (56) corresponding to any order of 
perturbation theory for The same is valid for the denominator of Eq. (20), Eq. (24) for 
the mass, and the expression for the renormalization constant (23). If instead we had chosen 
$(P) to transform differently from Eq. (32), even the fact that the solution of Eq. (24) does 
not depend on P would be hidden. 



3. Degenerate case 

In the degenerate case there is more than one solution $ of the unperturbed bound 
state equation, Eq. (8), for a single unperturbed bound state mass M u . Assuming an r-fold 
degeneracy, we denote such wave function solutions as $j where j = 1,2,3, ... ,r. In this 
case the pole structure of the unperturbed Green function G U (P) is easily seen to be 

Gu{P) = l ^t^fK G b u {P). (33) 

u 

As for the non-degenerate case, we shall assume our wave functions to be covariant but 
not dependent on P 2 . The wave functions $j are, by the assumption of r-fold degeneracy, 
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linearly independent. Applying this fact to the pole structure of the identity G U G U 1 G U = G u , 
we obtain the normalization condition for these wave functions: 



1 OP 2 3 



= 5ij. 



(34) 



P = Pu 



Eq. (33) can be written exactly as Eq. (7) with <£> now defined to be a row matrix whose 



elements are the 



$ = ( $i $ 2 $3 • • • $r ) , 



(35) 



with <£> being the corresponding column matrix with elements $j. With this redefinition of $ 
and $, the above derivation for the nondegenerate case remains valid up until and including 
Eq. (22). In this way we obtain, for the degenerate case, that 

G(P) = iil>(P)A-\P)$(P) + G\P), (36) 

where ip and ip are row and column matrices defined by elements 

UP) = [l + G 6 (P)^(P)] $,(P) and ^j(P) = $j(P) [l + ^(P)G 6 (P)] , (37) 

respectively, and A is an r x r matrix whose elements are 

Ay(P) = (P 2 - M 2 u )5 l3 - ih{P) [MP) + Ki{P)&{P)K x {P)] $,(P). (38) 

We are interested in the masses M for which the Green function G(P) of Eq. (36) 
develops a bound state or resonance pole. This will happen when the determinant of matrix 
A(P) becomes zero. This, in turn, can be determined by finding the matrix S(P) which 
diagonalizes A(P). With S(P) determined, we have that 



D{P) = S- 1 (P)A(P)S(P) 



and 



where 



and 



/P>i(P) o ••• \ 

D 2 (P) ••• 

D 3 (P) ■■■ 

V D r (P) J 

G(P) = i^j s {P)D- 1 {P)^ s {P) + G\P) 



D ii (P) = {P i -Mt)8 ij -i& i {P) K 1 (P) + K 1 (P)C?'(P)K 1 (P) *J(P), 



(39) 



(40) 



(41) 



U(p) = 4>(p)s(p), U(P) = s-\p)^(p), 



(42) 
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with similar definitions holding for $ S (P) and & S (P). Since detD(P) = T7j Dj(P) = 0, the 
Green function G(P) will have poles at P 2 = M 2 , j = 1, 2, 3, . . . , r, where Mj is the solution 
of the equation 

M] = Ml + i$f (P,-) [^(P,) + ^(^(P,)!^)] $f (P,), (43) 

Pj being any momentum satisfying Pj = Mj, and the functions $f(P) and $f(P) being 
the j'th elements of $ 5 (P) and $ 5 (P), respectively. 

Taking into account the diagonal nature of D(P), Eq. (40) can be written as 

G(P) = iY,^{P)D]\P)^{P) + G\P). (44) 

j 

Thus, assuming that the perturbed bound state mass Mj is itself nondegenerate, we can find 
its corresponding wave function tyj as in the nondegenerate case above: tyj = ^JZjijjj (Pj) , 
where 

Z ' J " 1 - i {If (P 3 ) \MP~) + K^G^K^)] if (P,)}" (45) 

Thus, in the degenerate case of the unperturbed theory, the properly normalized wave func- 
tions corresponding to the (nondegenerate) bound state mass Mj of the full perturbation 
theory, are 



^ = ^ [1 + G^Pj^Pj)] ^(Pj), (46) 
^ = $f (P,) [1 + K 1 (Pj)G b (Pj)}. (47) 



^. Comments 

The main results of this subsection are the expressions for M 2 and \& given in the 
nondegenerate case by Eq. (24) and Eq. (25), and in the degenerate case by Eq. (43) and 
Eq. (46), respectively. Not only are these expressions exact and compact, but they can also 
be easily used to write down the explicit perturbation series for these quantities. For this 
purpose it is most convenient to treat all functions of P as functions of P 2 and the unit four- 
vector n — P/ VP 2 , and at the same time to use the covariant form for the unperturbed 
wave function given by Eq. (31), as then $ will not depend on P 2 . For example, in the 
nondegenerate case, to generate the perturbation series for M 2 we use Eq. (12) to write 
Eq. (24) as an infinite series 



M 2 = M 2 + i$ 



K x + K X G U K\ + K\G U K\G U K\ + K\G U K\G U K\G U K\ + . . . $ (48) 



where a tilde over K\ or G b u indicates that this quantity is evaluated at P 2 = M 2 . By 
making Taylor series expansions 
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K 1 = K 1 + 6K[ + -K'{ + 



G b u = G b u + 5G b u + -G b a + ... 



where 



6 = M 2 - M 2 



(49) 
(50) 

(51) 



and each term without a tilde is evaluated at P 2 = M 2 , we can immediately write M 2 as a 
perturbation series with respect to orders of K 1 = Ki(M u ): 



M 2 = M 2 + 61 + 62 + 63 + ... 



where 



6 1 = iQK^ 



6 2 = i$ 

63 = i<5 
etc. 



6^ + KiG b u K^ $ 

6 2 K[ + 5 -±K'l + 6 1 (K&lK^)' + K\G\K\G h u K\ 



(52) 

(53) 
(54) 

(55) 



Similarly, the wave function of Eq. (25) can be written as a perturbation series in orders of 



= (1 + If! + 7] 2 + T]3 + . . .) $ 



(56) 



where 



Vl = |A X + G b u K 1 

V2 



(57) 

±A 2 + |A? + 5 1 (GlK,)' + ^G\K X (58) 

m = iAs + fAxAa + iA? + (fc + <W) + §<5? (G^i)" + ^G h u K x (59) 

etc. 

where Aj is derived from Si by putting an extra derivative on each K\ and G h u \ that is, 

Ai = i$K[$ (60) 



A 2 = i$ 

A 3 = i<J 
etc. 



<W + f + 6 X (K 1 (%K i y + (K^K^KJ 



(61) 
(62) 
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A similar procedure can be used to generate the perturbation series for the degenerate case. 

It is worth noting that the perturbative corrections to the bound state wave function, as 
derived here, are particularly important to take into account when calculating corrections 
to vertices (electromagnetic, axial, etc.) within constituent models. It is only by taking into 
account the appropriate order of wave function perturbation exactly, will symmetry proper- 
ties, like for example gauge invariance, be preserved at each order in the vertex correction 
- for a concrete example, see Ref. [10] where Eq. (57) was used to determine the full lowest 
order pionic correction to the nucleon vertex function in the NJL model. 

It is also worth pointing out that in the case where the perturbation K 1 is too large for 
a perturbative treatment, our expressions of Eq. (24), Eq. (25) Eq. (43), and Eq. (46) may 
still be useful for performing practical nonperturbative calculations of M 2 and Indeed, 
in both the degenerate and nondegenerate cases, the main calculational effort would be in 
solving Eq. (12) for the "background" Green function G b . Yet this is an especially simple 
equation, of standard Lippmann-Schwinger form, where G b has no pole at P 2 = M 2 and 
G\ has no pole at P 2 = M 2 (since they have been subtracted), and where there are no 
singularities in the integration over momenta. Even in the unlikely event that G b u happens 
to have an unsubtracted pole close to P 2 = M 2 , this case can be easily handled numerically. 
Finally, it is useful to note that G b u has already been constructed for the important case of 
the nonrelativistic Coulomb problem by Schwinger [11] - a result that can be easily adapted 
to the relativistic Coulomb case [2]. 

III. DISCUSSION AND SUMMARY 

In this work we have presented a general formulation of perturbation theory applicable 
to bound states and resonances where the bound state equations involve kernels and inverse 
free Green functions that have an arbitrary energy dependence. Our formulation is thus 
directly applicable to the important case of relativistic quantum field theory. One can 
consider our results as extending the well-known time-independent perturbation theory of 
quantum mechanics to the case where the kernels are energy-dependent and where the inverse 
propagators are non-linear in the energy. 

In particular, we have derived expressions for the bound state (or resonance) mass M 
and wave function ^ of a system whose interaction kernel K consists of a part K for which 
the corresponding Green function G u is known, and a part K 1 which plays the role of a 
perturbation. Our results for M and \1/ are contained in Eq. (24) and Eq. (25) for the 
nondegenerate case, and in Eq. (43) and Eq. (46) for the degenerate case, and have the 
feature that they are exact, with the perturbation K\ taken into account to all orders. The 
key element in these expressions is the Green function G b which needs to be found by solving 
Eq. (12). For sufficiently small Ki, Eq. (12) can be solved simply by iteration, in this way 
generating a perturbation expansion in K\ that is the analogue of the time-independent 
perturbation theory of quantum mechanics. On the other hand, if K 1 is not small enough 
to generate a convergent perturbation series, Eq. (12) could still be solved by standard 
numerical techniques for integral equations. 

As far as we know, our formulation of the perturbation theory problem is new. How- 
ever, there are a few alternative formulations available in the literature, all presented for 
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the particular case of relativistic quantum field theory. The first of these is a method where 
the perturbation series for M 2 and \I/ are expressed in terms of contour integrals. Originally 
developed by Kato [9] and described in Messiah's standard text [12] for the case of quantum 
mechanics, the contour method was extended to the covariant case by Lepage [5] and used, 
for example, by Murato [13]. Another method, due to Bodwin and Yennie [6], is closest in 
spirit to our approach, but does not have the feature of having closed expressions for the 
perturbed mass and wave function. A third approach is the recent formulation of Ivanov et 
al. [2] whose perturbative expansion is expressed in terms of a certain "relativistic general- 
ization of a projection operator" . In this approach the second derivative of the inverse free 
propagator, 8 2 Gq 1 /dE 2 , looks very much like a genuine and necessary relativistic feature, 
yet it does not appear in our formulation at all and is thus just an artifact of the particu- 
lar derivation used. Similarly, the expression for the lowest-order wave function correction 
derived directly from Eq. (9) of Ref. [2] contains four terms against our only one. 

In each of the above three alternative approaches, perturbative corrections to the bound 
state wave function were derived only for the special case where the bound state is at rest. 
Thus, in order to describe scattering process where the bound state has non-zero total 
momentum, such wave function corrections need to be modified by the appropriate Lorentz 
boost (that itself depends on the order of perturbation being considered). By contrast, our 
approach has enabled us to write expressions for the bound state wave function corrections 
that are Lorentz covariant at each order of the perturbation, thus avoiding the step of 
boosting from the rest frame. Although all perturbation expansions must mathematically 
be identical, it is evident that the expressions provided by our Eq. (24), Eq. (25), Eq. (43), 
and Eq. (46) are the simplest both practically and conceptually. 
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